Two-dimensional Non-Hermitian Derealization Transition as a Probe for the 

Localization Length 
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When one applies a type of non-Hermitian effect, constant imaginary vector potential, to disor- 
dered systems, derealization is induced even in two or lower dimension. By using the non-Hermitian 
induced transition as a probe, We propose a new procedure of estimating localization in arbitrary- 
dimensional systems. By examining numerically the two-dimensional non-Hermitian tight-biding 
model with onsite disorder, it is shown that the failure of absorbing the non-Hermitian effect, 
namely the breakdown of the imaginary gauge transformation, characterize the inverse localization 
length near the band center. 
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In non-interacting disordered systems with time- 
reversal symmetry, all the states are believed to be lo- 
calized in two or lower dimensions. Often one encounters 
situations where non-Hermitian effect emerges effectively 
through coupling with the reservoir, interaction between 
particles, or non-equilibrium relaxation processes. This 
type of non-Hermitian effect is taken into account as a fi- 
nite inelastic scattering time, hence following the scaling 
theory idea, localization will be truncated halfway only 
to exhibit a gradual crossover into metallic behavior. Re- 
cently, in a vortex pinning problem in a bulk supercon- 
ductor by columnar defects, it has been stressed that the 
non-Hermitian term coupled with the first-order differ- 
ential, i.e., the imaginary vector potential, gives a fun- 
damentally different effect on localization states in dis- 
ordered systemstl. The term appears to take account for 
the in-plane magnetic field component. Alppg with sub- 
sequent numerical and analytical studiesaQ, it has been 
recognized that sufficiently large constant imaginary vec- 
tor potential induces states with complex eigen ener- 
gies even in one dimension. The phenominon has later 
been understood as a non-Hermitian induced derealiza- 
tion transition not a crossover. The same class of non- 
Hermitian terms is ubiquitous in diverse problems: clas- 
sical (elastic, electromagnetic, or acoustic) wave propa- 
gation with absorptive mediumU, random walks in ran- 
dom media with drifting forcestl, or_Burgers turbulence 
through a Cole-Hopf transformationLj. 

While non-Hermitian effect is usually hard to control, 
the strength of the imaginary vector potential is rela- 
tively easy to regulate not only in numerical studies but 
also in experiments. For instance, one changes the tilt an- 
gle of the magnetic field in the vortex pinning problem, 
or tunes an incident frequency and change the absorp- 
tion rate in electromagnetic wave propagation. It has 
been suggested that the non-Hermitian induced dereal- 
ization process can be detected by the energy spectrum. 
So with its drastic consequence on localization, the phe- 
nomenon provides a great opportunity to investigate the 



localization property without resorting to any transport 
phenomenon, i.e., neither attaching leads nor coupling 
with the reservoir. 

In this paper, pursuing the above possibility, we pro- 
pose a new procedure (Eq. (||) below) appropriate to 
extract the inverse localization length out of the non- 
Hermitian dependence of the energy spectrum. Whereas 
the connection between the average density of states 
(DOS) and the localization length has been success- 
fully managed in ..one dimension, it is not the case in 
higher dimensionsa (See the argument below). Our pro- 
posed procedure is principally applicable in arbitrary- 
dimensional systems. By using the two-dimensional 
(2D) non-Hermitian model with time-reversal symme- 
try, we will confirm our proposal numerically. The non- 
Hermitian derealization in 2D systems is far from com- 
plete understanding, and it is stUL-at issue whether the 
transition is continuous or abruptBtx On this regard, we 
also briefly mention our point of view based on numerical 
results. 

The Hamiltonian considered for numerical calculations 
is the d-dimensional non-Hermitian tight-binding model, 

ffC0 = -f£ expiry) \i){j\ + J2 V ^^> C 1 ) 

where h is the directed hopping, i.e., a lattice analog of 
the imaginary vector potential. {V{\ is a random on-site 
potential chosen from a box distribution over (-W, W). 
\i) denotes the state bound at r^, and the lattice sum- 
mation (i,j) is restricted over nearest- neighboring pairs 
of sites. We work on the two-dimensional case, where 
we particularly choose h as h x = h y = h to avoid un- 
necessary troubles from lattice periodicity. In the rest 
of the paper, we adopt the convention of the unit lattice 
constant and t = 2 as well as % = 1. Using an exact di- 
agonalization technique, we have calculated a complete 
set of eigen values and eigen wavefunctions on a two- 
dimensional 20 x 20 lattice. (Eigen functions are needed 
later to justify our interpretation.) Ensemble average is 
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taken over 40 samples. 

The possibility to detect localization by using the non- 
Hermitian dependence of the spectrum has already been 
noticed in Refill. The reasoning is as follows. Suppose an 
eigen state (j)h=o(f) be localized exponentially around ro 
with the localization length £, that is, 0o ~ e~l r ~ ro l/'». 
When one turns on h gradually, its eigen energy does 
not change in value, while the right eigen wave func- 
tion progresses following the imaginary gauge transfor- 
mation 4>h{f) ~ e h '(' r °^4>a{r) up to a small correction 
by the boundary condition. Such 'absorption' of the non- 
Hermitian influence on the spectrum can take effect only 
when h is smaller than the inverse localization length. 
Beyond that, the state undergoes such a drastic change 
that it transforms into a complex energy (delocalized) 
state. Putting the other way round, one uses the break- 
down of the imaginary gauge transformation to charac- 
terize and give a practical estimate of the inverse local- 
ization length. 

A naive way to estimate the breakdown of the imag- 
inary gauge transformation is to examine the average 
DOS defined on the complex energy plane, ph(E,£) = 
(6(E — ReE n )S(£ — lmE n )). When all the states are lo- 
calized at sufficiently small h, energies are all real and 
Ph(E,£) oc S(£). With increasing h, the support of 
Ph{E,£) extends into a £ ^ region in the presence of 
complex energy states. By defining ho(E) by the point 
having a finite density at £ ^ for each E one sees 
ho(E) indicate the breakdown of the imaginary gauge 
transformation. Accordingly, it is suggested that the in- 
verse localizatiottlength may be estimated by the formula 

i/t(E) = h (E)m. 

While the above argument is conceptually simple and 
attractive, the estimate of the localization length by us- 
ing the average DOS has been successful only in one 
dimension, whose success is traced to a special struc- 
ture of the energy spectrum. In one dimension, irrespec- 
tive of the value of h, the support of ph{E,£) becomes 
one-dimensional on the complex energy plane, i.e., con- 
sisting of two line segments along the real axis plus a 
'bubble' part, in-between corresponding to complex en- 
ergy statesElia. The shape of the bubble turns out quite 
insensitive to disorder configurations. Because of it, one 
can readily read off the value of h$ (E) from the endpoints 
of the arcs of the bubble. 

By contrast, the DOS support on higher dimensions 
is no longer one-dimensional but much more involved. 
A typical structure of the two-dimensional DOS for a 
sufficiently large h is shown in the inset of Fig. fil. As 
one sees, localized and delocalized states coexist unpre- 
dictably over the entire energy range, showing that DOS 
suffers large fluctuations. The system has no clear-cut 
mobility edge distinguishing between real-energy (local- 
ized) and complex-energy (delocalized) states.— .In fact, 
judging from existing numerically and analytically results, 
it is very likely that the average DOS p(E,£) always 
retains a finite contribution proportional to 5(£) hence 
real energy states persist in existing for sufficiently large 



h, even near the band center. The feature of ph(E,£) 
in higher dimensions makes it highly dubious to seek 
the breakdown of the imaginary gauge transformation 
at a fixed E. Moreover, numerical calculations on two- 
dimensional systems suggest that ho(E) will estimate an 
inverse localization length at more than one order of mag- 
nitude larger than that obtained by any other method 
such as the recursion methocO. 

Resolving problems above to answer whether the non- 
Hermitian effect characterizes the localization length or 
not, we propose an alternative criterion to extract the lo- 
calization length. Our solution is to find the breakdown 
of the imaginary gauge transformation in a more faithful 
and careful way: First we identify the smallest directed 
hopping strength h c where any of the eigen energies gets 
complex in a certain configuration of the random poten- 
tial (see Fig. [I]) , and next we take the ensemble .average 
of h c to obtain the inverse localization length byEJ 

(h c ) = I/£(0). (2) 

The estimated localization length corresponds to that 
around the band center, since states near the band center 
tend to be delocalized first. Eq. (Q) and the procedure 
leading to it is our main result in the paper. The advan- 
tage of Eq. @ is clear. In any localized system, we can 
assign an unambiguous value of h c for each configuration, 
and perform the ensemble average of h c . It should be 
stressed that although there is no discrepancy between 
(h c ) and /io(0) for one dimension, (h c ) becomes much 
smaller than /io(0) in higher dimensions, since the lat- 
ter cannot detect a weakly-broken situation as in Fig. [|. 
Unlike ho (E) (if at all correctly defined in higher dimen- 
sions), we cannot obtain (h c ) from the average DOS. 

It is worth discussing a little more the definition of 
h c and finite volume effect. Although we have initially 
intended to assign h c as the point where the imagi- 
nary gauge transformation fails to absorb all the non- 
Hermitian effect, we should keep in mind that such mech- 
anism works only approximately in a finite system. Fig- 
ure |^ represents a typical evolution of the first pair of 
levels getting complex. It shows that complex energy 
states appear through merging two adjacent levels, but 
well before it, the pair of levels shift gradually without 
gaining any imaginary part. Strictly speaking, the effect 
of h remains partially for any small h, hence the gauge 
transformation is already broken down for infinitesimal 
h. Nevertheless, we assign the value of h c by a merging 
point. The shift of energies before h c arises from finite 
size effect, the consistency with the periodic boundary 
condition in a finite system. 

To justify the statement above, Fig. || illustrates the 
spatial dependence of the right eigen wave function <fi(r) 
of the corresponding levels of Fig. [j]. When we look at 
the amplitude next to a localization center, it evolves 
precisely in the form predicted by the imaginary gauge 
transformation (the solid line in the figure). However, 
such dependence fails appreciably at a point sufficiently 
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far away from a localization center (r = 9 on a 20 x 20 
lattice). The deviation is responsible for the level shift 
observed in h < h c , and has nothing to do with any pre- 
cursor to derealization. Simultaneously, Fig. || reveals 
clearly that it is almost impossible, at least numerically, 
to detect the point h c by examining the spatial depen- 
dence of the wavefunctipn amplitude, in contrast to the 
one-dimensional systemEj. The overall shape is barely 
changed through the critical value h c since the contribu- 
tion of the imaginary part is too small to give a visible 
deviation (see the inset). It is remarked that on a few 
occasions some abrupt movement of the localized center 
is observed before reaching h c without any finite imagi- 
nary part. We still find, however, the spatial dependence 
following the gauge transformation close to a localiza- 
tion center both before and after it, so the phenomenon 
should be attributed as well to an artifact of finite system 
sizetS 

To obtain h c from numerical data, we need to distin- 
guish between energies with vanishing and non- vanishing 
imaginary parts. Such type of estimate is usually hard 
in numerics, but in our problem, we can unambiguously 
identify it, because complex energy states always emerge 
by pairs (see Fig j^) . 

Following the procedure leading to Eq. (^) , we find h c 
for a given configuration of the random potential, and 
take the ensemble average to obtain the inverse localiza- 
tion length. Our numerical results are compared with the 
ones obtained from the. recursion method combined with 
finite scaling analysisliJ, which is known among the most 
reliable estimate of the localization length. Since our nu- 
merical parameters reside in a strongly localized region, 
we take the Lyapunov exponent directly as the inverse 
localization length with little finite size correction. 

Figure |] shows the ensemble-averaged value (h c ) and 
its variance from a 20 x 20 lattices, compared with the 
inverse localization length for aooxoo lattice denoted by 
the solid line (for— L6 x oo lattice by the dashed line) ex- 
cerpted from Rcfo. The agreement between (h c ) and the 
inverse localization is clear, and the variance 5h c is rela- 
tively small, considering the fact we work on a compar- 
atively small system size and ensembles. This shows the 
effectiveness of the present approach. Our data points 
show a systematic drift, though relatively small, from the 
oo x oo line to the 16 x oo line in approaching h c — > 1. The 
tendency is consistent because finite size effect becomes 
increasingly prominent in this limit. It is emphasized 
that the agreement is substantial and the ability to pre- 
dict the localization length quantitatively is the virtue of 
Eq. (||), while other numerical objects such as the par- 
ticipation ratio are often lacking. 

As for the nature of the two-dimensional non- 
Hermitian derealization transition, our view based on 
the present numerical results is as follows. By attribut- 
ing the level shift in the region h < h c to finite size effect, 
it is natural to infer that the transition is not continuous 
but abrupt in the limit of infinite size at least for a given 
configuration. By assuming further the self-averaging 



property, one is tempted to reach the conclusion that the 
two-dimensional non-Hermitian derealization transition 
is abrupt for bulk systems. It should be recalled, how- 
ever, that in our particular choice of averaging procedure, 
we hardly respect the additional symmetry present at the 
band center. So we believe as a possible scenario, that the 
transition away from the band center is discontinuous but 
continuous at the band center, which is consistent both, 
with the nonperturbative analysis by the supermatrix£L2l 
and with the symmetry argument applicable to the band 
centena. However, the present results are not sufficient 
to conclude this clearly, and it requires further studies in 
the future to draw a decisive conclusion on the issue. 

In conclusion, we have shown for the first time how the 
two-dimensional localization can be probed by using the 
non-Hermitian transition induced by the imaginary vec- 
tor potential, proposing a practical procedure to make an 
estimation of the inverse localization length. By checking 
the level shifting through the transition, it has been un- 
derstood clearly that the two-dimensional non-Hermitian 
derealization transition occurs through merging adja- 
cent levels. It is true at present that the evaluation of the 
localization length by using the non-Hermitian effect has 
not yet fulfilled such accuracy as the recursion method 
enjoys, and it needs much more labor or fine tuning to 
find h c - However, it surely serves as a useful alterna- 
tive that one may be able to apply even in experimental 
settings. Furthermore, as the present method need not 
increase the system size, it may be advantageous in nu- 
merical studies on more complicated systems such as the 
disordered interacting systems. 

We acknowledge the support in the initial stage of 
the research from the Grant-in- Aid for Scientific research 
No. 11216204 by Japanese Ministry of Education. 
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FIG. 1. An typical example of the energy spectrum for a 
given configuration of random site energies (only the upper 
half plane shown here). Random on-site potential is chosen 
from a box distribution [—10, 10]. The value of h is adjusted 
to h — 0.21, where only one pair of states (indicated by the 
arrow) get complex. Inset: A typical structure of the DOS for 
a larger h (h — 0.60). Both localized and delocalized states 
are mixed up in the same energy range. 




FIG. 2. An evolution of the first pair of levels turning 
complex values, by increasing h. In this configuration, we as- 
sign the critical value h c as h c = 0.21. The lmE = plane is 
drawn as a shaded plane for the eye guide. 




FIG. 3. The validity of the imaginary gauge transformation 
in the right wave function. r _1 log |Re<l?(ro +r)| is plotted as 
a function of the imaginary vector potential h where ro is the 
localization center. The value is normalized by that at h — 0. 
Squares refer to the values at r = 1, and open triangles to 
those at r — 9. A solid line is what is expected by the imagi- 
nary gauge transformation. Inset: |Im$(ro + r)\ for r — 1 as 
a function of h. A jump at h c is clearly observed. 




FIG. 4. (h c ) versus the disorder strength W to confirm 
Eq. (Q) . Data points (o) are obtained from 40 ensembles of a 
20 x 20 lattice. The solid line is the inverse localization length 
for aooxoo latticer-and the dashed line for a 16 x oo one. 
Excerpted from Ref.tJ. 
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